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ABSTRACT 

We investigate the practical implementation of Taylor's (2002) 3-dimensional grav- 
itational potential reconstruction method using weak gravitational lensing, together 
with the requisite reconstruction of the lensing potential. This methodology calculates 
the 3-D gravitational potential given a knowledge of shear estimates and redshifts for a 
set of galaxies. We analytically estimate the noise expected in the reconstructed grav- 
itational field taking into account the uncertainties associated with a finite survey, 
photometric redshift uncertainty redshift-space distortions, and multiple scattering 
events. In order to implement this approach for future data analysis, we simulate the 
lensing distortion fields due to various mass distributions. We create catalogues of 
galaxies sampling this distortion in three dimensions, with realistic spatial distribu- 
tion and intrinsic ellipticity for both ground-based and space-based surveys. Using the 
resulting catalogues of galaxy position and shear, we demonstrate that it is possible 
to reconstruct the lensing and gravitational potentials with our method. For example, 
we demonstrate that a typical ground-based shear survey with redshift limit 2 = 1 
and photometric redshifts with error Az = 0.05 is directly able to measure the 3-D 
gravitational potential for mass concentrations ^ 10 14 M© between 0.1 z ^ 0.5, and 
can statistically measure the potential at much lower mass limits. The intrinsic ellip- 
ticity of objects is found to be a serious source of noise for the gravitational potential, 
which can be overcome by Wiener filtering or examining the potential statistically 
over many fields. We examine the use of the 3-D lensing potential to measure mass 
and position of clusters in 3-D, and to detect clusters behind clusters. 

Key words: Gravitation; Gravitational Lensing; Cosmology: Observations, Dark 
Matter, Large-Scale Structure of Universe. 



1 INTRODUCTION 

Gravitational lensing affords us a direct method to probe 
the distribution of matter in the universe, irrespective of its 
state or nature. This deflection of light by the gravitational 
potential of matter along its path can be observed as a local 
alteration of number counts of background galaxies (magni- 
fication), or a distortion of their shape (shear). It is the latter 
phenomenon which we will concern ourselves with here; we 
will further restrict ourselves to the case where this distor- 
tion is weak (^ 10% change in the ellipticity of the object). 
Despite the weakness of the effect, and the intrinsic, nearly 
randomly orientated ellipticity of background galaxies, we 
can measure the weak shear by averaging the ellipticity or 
shear estimates of very many galaxies. 

It has long been recognised that weak gravitational lens- 
ing is a valuable tool for examining the two dimensional pro- 
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jected matter distribution, and can consequently provide im- 
portant information regarding large-scale structure (see e.g. 
Bartelmann & Schneider 2000, Bernardeau et al 1997, Jain 
& Seljak 1997, Kaiser 1998). In particular, the sensitivity 
of lensing to all the matter present, including the dominant 
dark matter, ensures that the weak lensing effect is an ex- 
cellent probe for determining the quantity and distribution 
of matter. 

Weak lensing studies for a wide range of galaxy clusters 
have been carried out, allowing precision measurements of 
the clusters' masses and mass distributions (see e.g. Tyson et 
al 1990, Kaiser & Squires 1993, Bonnet et al 1994, Squires 
et al 1996, Hoekstra et al 1998, Luppino & Kaiser 1997, 
Gray et al 2002). Moving to larger scales, the shear due 
to large-scale structure has been accurately measured by 
several groups (see e.g. van Waerbeke et al 2001, Hoekstra 
et al 2002, Bacon et al 2002, Refregier et al 2002, Brown et 
al 2002). 

Redshift information has already been used in weak 
lensing studies, e.g. to determine the median redshifts of 
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the lens and background populations; most analyses then 
project the lensing information into a 2-dimensional pro- 
jected mass distribution. Wittman et al (2001, 2002) demon- 
strate the utility of using 3-D shear information by inferring 
the redshift of a cluster using shear and photometric redshift 
information for the galaxies in their sample. The importance 
of including redshift information to remove intrinsic galaxy 
alignments from shear studies has also been discussed (Hey- 
mans & Heavens 2002, King & Schneider 2002a,b). Lens 
tomography has been studied as a valuable means of intro- 
ducing redshift information into shear power spectra (e.g. 
Seljak 1998, Hu 1999, 2002, Huterer 2002, King & Schnei- 
der 2002b), but only recently has the full reconstruction of 
the 3-D Dark Matter distribution from lensing been consid- 
ered (Taylor 2002, Hu & Keaton 2002). 

In this paper, we seek to discuss a practical implemen- 
tation for reconstruction of the 3-D lensing and gravitational 
potentials from weak lensing measurements and redshift in- 
formation (whether photometric or spectroscopic). We will 
use the reconstruction procedure of Taylor (2002) , which al- 
lows us to calculate the entire 3-D gravitational potential if 
we have a knowledge of shear estimates and redshifts for a 
set of galaxies. This procedure is explained in Section 2. 

In Section 3 we discuss sources of uncertainty for our re- 
construction, examining analytically the effect of shot noise 
due to galaxy ellipticity, the effects induced by a finite sur- 
vey, photometric redshift errors, redshift-space distortions, 
and multiple scatterings of light rays. 

We aim to test the reconstruction method using simu- 
lations of realistic weak lensing data. In practice, for a par- 
ticular volume of space containing a mass distribution, the 
data which would result from a real survey would be a set of 
galaxy ellipticities (in the weak lensing regime, these will be 
dominated by intrinsic ellipticity with a small gravitational 
lensing perturbation) and redshifts. The ellipticities will typ- 
ically be defined by the galaxies' quadrupole moments (e.g. 
Kaiser et al 1995, Rhodes et al 2000) or estimated from a de- 
composition of galaxy shape into eigenfunctions (Refregier 
& Bacon 2002, Bernstein & Jarvis 2002). In order to pro- 
vide a realistic catalogue of these data, we simulate mass 
distributions, calculate the expected lensing distortion in 3- 
D, and create a set of galaxies probing this distortion with 
a realistic redshift distribution and intrinsic ellipticity. We 
describe this procedure in detail in Section 4.1. 

Given such a catalogue, we attempt to reconstruct the 
lensing potential and gravitational potential (Section 4.2). 
We use a generalised 3-D Kaiser-Squires (1993) inversion 
together with Taylor's (2002) formalism to obtain the 3- 
dimensional gravitational potential distribution. 

In Section 5 we demonstrate the effectiveness of our im- 
plementation in reconstructing lensing and gravitational po- 
tentials. We examine the level of noise in our simulations in 
Section 6, including the Poisson noise associated with having 
only a finite number of galaxies probing the lensing distribu- 
tion, and the noise from galaxies' intrinsic ellipticities. We 
go on to examine the uncertainties caused by photometric 
or spectroscopic redshift errors. 

In Section 7 we study the utility of reconstructing only 
the 3-D lensing potential without continuing to the gravita- 
tional potential; this provides a useful method for detecting 
mass concentrations and measuring their mass and 3-D po- 
sition. Finally, we summarise our results in Section 8. 



2 3-D GRAVITATIONAL POTENTIAL 
RECONSTRUCTION 

In this section we will summarise the results of Taylor's 
(2002) approach to reconstructing the 3-D gravitational po- 
tential, using weak lensing measurements together with red- 
shifts for all galaxies in the lensing catalogue. 

We begin by noting that the Newtonian gravitational 
potential $ can be related to the density of matter p by 
Poisson's equation, 



V 2 $ = 47vGp m 8a 2 — -\ H 2 £l m a 1 S, 



(1) 



where we have introduced the cosmological scale factor a, 
the density contrast S = (p — p)/p, the Hubble length 
Xh = l/Ho ~ 3000 ft -1 Mpc, and the present-day mass- 
density parameter f2 m . 

We can study the impact of this gravitational poten- 
tial on image distortion due to gravitational lensing, by in- 
troducing the lensing potential <f>. This is a measure of the 
distortion which is related to the observable shear matrix by 



Hi 



(2) 



where di = r(5ij — fifj)Vj = r(Vi — fid r ) is a dimensionless, 
transverse differential operator, and d 2 = did 1 is the trans- 
verse Laplacian. The indices take the values (1, 2), and 
we have assumed a flat sky. 

The lensing potential is also observable via the lens con- 
vergence field 

k = \d 2 <t>. (3) 

The convergence field is related to the shear field by the 
differential relation, first used by Kaiser & Squires (1993); 

k = d~ 2 didjiij, (4) 

where d~ 2 is the inverse 2-D Laplacian operator on a flat 
sky, defined by 



d~ 2 ee — / d 2 9 ln|0-0'| 
2tt J 



(5) 



Note that transverse positions 6 in this equation will be 
quoted in units of radians, so that the lensing quantities 
are dimensionless. However, such angular positions can be 
arbitrarily scaled, leading to simple scalings on the lensing 
quantities which we will quote later. 

A useful quantity for tracing noise and systematics in 
gravitational lensing is the divergence-free field, /3, defined 
by 



P = d 2 e?d J d„"f ij , 



(6) 



where e? 



[ f is generated purely by the 



-1 

1 

lensing potential, f3 vanishes. But if there are non-potential 
sources, due to noise, systematics or intrinsic alignments, 
then p will be non-zero. In addition, P terms can arise from 
finite fields, due to mode-mixing of shear fields (e.g. Bunn 
2002). We discuss this further in Section 3.2. 

We can relate the lensing potential and the gravitational 
potential by 



</>(r) 



2 / dr 



rr' 



$(r') 



(7) 



© 2002 RAS, MNRAS 000, 1-21 



Mapping the 3-D Dark Matter Potential 3 



in a spatially flat universe, with comoving distance r. This 
equation assumes the Born approximation, in which the 
path of integration is unperturbed by the lens. 

Note that the lensing potential is really a 3-D quan- 
tity, although usually it is regarded as a 2-D variable in the 
absence of redshift information. In this case the usual prac- 
tice is to average over the redshift distribution of the source 
galaxies (e.g. Bartelmann & Schneider 2001). It is in this 2-D 
approximation that the depth information is lost in lensing. 

Given that 4>{r) is really a 3-D variable, we can readily, 
and exactly, invert equation (7) and recover the full 3-D 
Newtonian potential (Taylor 2002); 

*(»") = \d r r 2 d r <t>{r) (8) 

where d r = f.V is the radial derivative. Any lensing field 
from real data will contain significant noise, and will thus 
require smoothing if we are to perform this differentiation. 
Interestingly it turns out that the lensing potential obeys 
a second-order differential equation, which is r 2 times the 
radial part of the 3-D Laplacian. This appears to be just a 
coincidence given that the lensing kernel in equation (7) is 
solely due to the geometric properties of the lens. 

In order to reconstruct the gravitational potential, it is 
necessary to find the lensing potential from the shear. This 
can be achieved by using the Kaiser-Squires (1993) relation, 
generalised to 3-D: 



(r) =2d~ 4 d i d j 7y (r), 



(9) 



where <j> is an estimate of <j>- As the variance of the shear 
field is formally infinite (Kaiser & Squires 1993), this distri- 
bution is usually binned and/or smoothed in the transverse 
direction before calculating the lensing potential. 

The above solution allows us to estimate the lensing po- 
tential only up to an arbitrary function of the radial distance 



Similar relationships between the potential and the 3-D 
lens convergence can also be written down: 

$ = d r r 2 d r (d- 2 n - (<9~ 2 k», (13) 

while the relationship between the matter density field and 
the 3-D convergence is 



V d r r d r (d k — (d k)). 



(14) 



With these sets of equations, the 3-D lensing convergence, 3- 
D lensing potential, 3-D Newtonian potential and 3-D mat- 
ter density fields can all be generated from combined shear 
and redshift information. 



3 THE UNCERTAINTY IN 3-D LENSING 

3.1 Shot-noise uncertainty in lensing fields 

Having written down the basic equations for the 3-D analysis 
of gravitational lensing data, we now consider the various 
contributions to the uncertainty in a reconstruction of these 
fields. 



3.1.1 The convergence field 

The covariance on a reconstructed, continuous convergence 
field due to shot-noise is generally given by 

7r 2 m 

{n{r)K{r'))sN = -j^-Soir - r'), 



n(r) 



(15) 



where 7 rms is the intrinsic dispersion of galaxy shear es- 
timates in one component (i.e. 71 or 72) due to the non- 
circularity of galaxies, and n(r) is the observed space density 
of galaxies in a survey. 

In the case of a discretised map this reduces to 



4>(r) = 4>{r) + tp{r), 



(10) 



where <f> is the true lensing potential, and ip(r) is a solution 
to 



{pidj - ^5^d 2 ^j tp = 0. 



(11) 



This arbitrary radial behaviour is due to the fact that the 
shear only defines the lensing potential up to a constant for 
each slice in depth. However, we must tame this behaviour 
if we wish to apply equation (8) to find the gravitational 
potential. 

Fortunately, there are several opportunities for remov- 
ing i\>. Firstly, since (<f>) = for each slice in depth from ho- 
mogeneity and isotropy conditions, we can simply subtract 
(<j>) from <j> for each radial slice, provided that our survey is 
large enough to discount errors in {<j>) . Thus 



$= ^d r r 2 d r ($- 



(12) 



Alternatively, we can note that (d r (j>) = for each radial 
slice, with a similar subtraction required. Given these pro- 
cedures, it is clear that ip can be removed for large surveys, 
where the necessary averages (<j>) or (d r (j>) will have small 
uncertainty; however for limited area surveys ip will not be 
estimated well, leading to uncertainties on the reconstruc- 
tion (see Section 3.1.2). 



/ \ 7rms ?K 
(KiKjjsN = — 0, 



N ■ 



where, for a constant 3-D number density of galaxies, 



N, 



plx = 3n8 pix r Ar/R 



(16) 



(17) 



is the 3-D pixel occupation number, n is the 3-D density of 
galaxies, 9 p i x is the pixel size, Ar is the width of the radial 
bins, r is the radial distance, and R is the chosen limiting 
distance for the survey. We will compare this amplitude and 
behaviour with our simulations in Section 6. 



3.1.2 The lensing potential field 

We now wish to describe the uncertainty expected for the 
lensing potential field. We can write the covariance of the 
lensing potential estimated from equation (3) as 

(ct>{r)ct>{r'))sN = 4d- 2 d'- 2 {n(r)n(r')) S N. (18) 

We describe the procedure used to evaluate this covariance 
in the Appendix, and here only quote the resulting variance 
of the lensing potential at the centre of the survey (0 = 0) 
as a function of survey size, 9, and radial position, r, in the 
flat-sky limit: 



<A^ = ))s „ = ^^ Mr -O. 



(19) 
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Note that the uncertainty on the potential difference in- 
creases with the survey area as a consequence of the 2-D 
flat-sky lensing "force" term, In \0 — 0'\, increasing with dis- 
tance. We will compare this behaviour with our simulations 
in Section 6. 

Finally we can cast this in a more convenient form for 
gravitational lensing in the discrete case: 



(A^ 2 (0)) sn = 7.6x10" 



ri2 



30/[l' 



r 2 Ar 



,(20) 



where we have assumed 7 rms = 0.2, n 2 = |n7? 3 is the 2-D 
total surface density of galaxies and R is the nominal depth 
of the survey. 

In general the uncertainty in the 3-D lensing potential 
will also depend on the pixel size as well as the total size of 
the survey. This behaviour for the uncertainty on the lens- 
ing potential is now due to the 2-D flat-sky lensing "force" 
term diverging at small separation. Due to the complexity 
of analysing this effect analytically we shall defer a more 
thorough treatment of studying pixelisation effects with our 
simulations to Section 6. 



3.1.4 The density field 

Finally the shot-noise uncertainty in the reconstructed den- 
sity field can be constructed from equation (1); 



{6(r)5(r')) 



2X 2 H a 



V 2 V' 2 ($(r)$(r')). 



(24) 



In the far-field approximation the Laplacian can be writ- 
ten V 2 = (d 2 + R~ 2 d 2 ). In general equation (24) must be 
calculated numerically, but for points along the centre of 
the survey the shot-noise contribution to the variance of the 
density field can be calculated analytically, reducing to 



( 2 Y (_a_\ 2 jLs /AhV 
V<>.. / nrf, V R J 



V97T 3 / 



HI 



6R 



175 

24 



9R 

HI 



(25) 



If we set a = 1/2, and 7 rms = 0.2, and again define 
n 2 — |ni? 3 as the surface galaxy density, the variance of 
the density field can be expressed in the form 



3.1.3 The Newtonian potential field 

We also wish to calculate the uncertainty on the 3-D New- 
tonian potential. We can approach this, from equation (8), 
by differentiating the lensing potential field leading in the 
far-field limit to 



<*(r)S(r')>sw = irV 2 d 2 d; 2 <0(r)^(r')}siv. 



(21) 



We again describe the detailed calculation of this quantity 
in the Appendix, where we arrive at an expression for the 
variance on the Newtonian potential: 



(4> 2 (r)) S iv = 



Trms 



64tt 



2tt r 3 n 



(22) 



where r\\ is the smoothing radius for a radial Gaussian 
smoothing of the field. Hence we find that the shot-noise 
uncertainty is a strong function of the radial smoothing, 
changing as the inverse fifth power of the smoothing radius. 
We will again examine this behaviour in our simulations in 
Section 6. 

Finally, we can again cast this in a more convenient 
form for lensing: 



(4> 2 (r)) S iv = l.lxl0- 16 



n 2 



\30/[V] 



(f) : 



Note that, as for the lensing potential, the error upon the 
gravitational potential grows with survey or cell size. 

The question of what survey area and radial smooth- 
ing to choose leaves us with an optimisation problem. From 
equation (22) it would appear that we can reduce the shot- 
noise by reducing the survey area, or increasing the radial 
smoothing. However the latter will reduce the resolution of 
the survey, and correspondingly lower the intrinsic signal, 
while the former will increase the noise effects induced by 
a finite survey area. We shall investigate further the effects 
arising from a finite survey in Section 3.2, and comment 
further on the problem of survey optimisation. 



(S 2 (r))sN 



2.9 x 10" 



Ah 

Hi 



03 



n 2 



[30/1'] 



HI 



R 3 



1 + 6.8 x 10" 



(0/l°)R 



HI 



(26) 



To illustrate this we find that the variance in the recon- 
structed density field for r\\ = 0.1, R = 1 and r = 1, where 
we express distances here in redshift for a Euclidean uni- 
verse, with ?i2 = 30 galaxies per sq. arcmin. and Q. m = 0.3 
is given by 



(S 2 (r)) SN = 0.003 



(f) 



1 + 0.007 ( — 
' 1° 



(27) 



which has a minimum at 9 = 3.2° of {S 2 } 1 ^ 2 = 0.02. Since 
we expect the amplitude of density perturbations on these 
scales, A « 100 /i -1 Mpc, to be smaller than this, ag w 0.01, 
we can expect that filtering (e.g. Wiener filtering, see Section 
3.3) will be required to extract a large-scale map of the 3-D 
density field, even at the scale which minimises noise. 



3.2 Uncertainty and mode-mixing due to finite 
fields 



.(23) 



3.2.1 Finite surveys 

As well as shot-noise arising from the discrete sampling of 
the shear field by the survey galaxies, there is an additional 
uncertainty in the reconstruction of the density field for fi- 
nite area surveys due to the reconstruction process being 
nonlocal (via the inverse Laplacian) . The nonlocal behaviour 
of density reconstruction over a finite survey area also gives 
rise to a mixing of modes. Here we try to quantify for the 
first time for lensing the effects of mode-mixing and recon- 
struction noise arising from finite survey areas. The following 
analysis will be applicable to either 3-D or 2-D lensing stud- 
ies, as the results will be true for either a series of redshift 
slices or an overall 2-dimensional projection. 
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We may calculate the effect of a finite shear sample by 
multiplying the shear field by an arbitrary window function; 



■y'ij (r) = W / (r)7 lJ (r). 



(28) 



The lens convergence field can be estimated via equation (4), 
but since here there is only a finite field to integrate over, the 
effect of the window function is to truncate the effects due 
to the distant shear field and induce mode-mixing between 
k and p. 

If we expand the shear in Fourier modes on the flat sky, 
the observed shear is a convolution of the intrinsic shear and 
the window function. For convenience we shall use continu- 
ous transforms, although for a finite sky a discrete Fourier 
transform with suitable boundary conditions is more prac- 
tical. A scalar quantity f(6,r) can be expanded in a 2-D 
Fourier series by 



Mr) 



Sl Mr)e* e . 



(2w) 



(29) 



Decomposing the observed shear matrix given by equation 
(28) into the Fourier decomposed k(£) and (3(£) fields using 
equations (4) and (6), and Fourier transforming again we 
find the relationships between the reconstructed and true k- 
and (5- modes are 

K 'W = / W(£-?)[K(£')cos2ip- f3(£')sm2ip], 
J (27T) 



(3'(£) 



W{£ - £')\f3{£') cos2<^ + k(£') sin 2^], (30) 



where cos tp — l.i . Hence we see that a finite survey will 
give rise to a spurious /3-field. Here we have suppressed the 
radial dependence for clarity, so these equations are directly 
applicable to reconstruction of the convergence field in 2-D 
lensing. 

If we assume negligible intrinsic (3 fields then the real- 
space convergence is given by 



K'(0) = 

where 
W t (p) = 



f dH K 



{£)W e (0)e 



-t 2 B 2 /2 



■W(l-l')cos2y e 



jL'.e 



(31) 



(32) 



is the total effect of the finite survey area, and we have 
assumed the shear field is Gaussian smoothed on a scale 8 S . 



3.2.2 Variance of the convergence field 

The variance measured in the finite-survey convergence field 
on a smoothing scale of 8 S is given by 



w 2 m 

where 



J^C t \W t (0)\ e - 



(X{£)Y{£)) = (2nfCf Y 5 D (£ - t) 



(33) 



(34) 



for the isotropic fields X(G) and Y{0). 

If we approximate the window function by a Gaussian 
of radius 8 with Fourier transform 



W{£) = 2n8 2 exp(-£ 2 /28 2 ), 



(35) 
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Figure 1. The finite-field lensing window function, Wt(0), for 
a Gaussian window of radius 9, calculated at the centre of the 
survey as a function of id. 



we can evaluate We(0) at 6 — 0, yielding 



Wi(0) = 1 + 



20 



,-t 2 e 2 /2 



1) 



£ 2 8 2 



(36) 



We plot this window function, giving the contribution of 
convergence modes to the observed convergence, in Figure 
1. The general effect of the window function is to act as 
a high-pass filter. At I ~ 1/8 all the convergence modes 
are destroyed, while at large 19 the window function tends 
to unity as all the modes on scales below the scale of the 
survey contribute to the variance. Interestingly at low 18 this 
function goes negative, indicating that modes larger than the 
survey area are heavily distorted. 

For the Gaussian window function, no /3-modes are gen- 
erated at the centre of the field due to the symmetry of the 
window. Hence the /3-modes that are generated by the finite 
window are distributed nonlocally over the survey area. 

Figure 2 shows the variance of the observed convergence 
field, k' , on a smoothing scale of 8 S — 0.05°, measured at 
the centre of a Gaussian survey as a function of the survey 
radius, 8 (solid line). We have assumed a convergence power 
spectrum, C^ KK , for a flat LCDM universe with Q m = 0.3, 
CIa = 0.7 and used the Peacock-Dodds transformation (Pea- 
cock & Dodds 1996) to map to the nonlinear regime. Below 
a survey radius of 8 = 0.1 deg., missing modes due to the 
finite window and mode-mixing result in a drop in the mea- 
sured variance. 



3.2.3 Uncertainty in the reconstructed convergence 

In addition to estimating the variance in the reconstructed 
convergence field (fi;'(0)}obs measured from within a finite 
survey, we can also predict the uncertainty («;'(<?))miss due to 
missing structure in the shear field beyond the survey area. 
As the total variance measured within the survey and the 
missing modes from beyond the survey must yield the total 
variance of the convergence field (/t'(0))totai, the uncertainty 
in a reconstruction due to missing structure is 
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[deg] 

Figure 2. Variance of reconstructed convergence field, k', Gaus- 
sian smoothed on a scale 8 S = 0.05°, measured at centre of Gaus- 
sian survey, as a function of survey radius, 6 (solid line). The 
variance in the reconstructed convergence due to shear structure 
beyond the survey boundary is also shown (dotted line). 

<K'(0)) miss = (/t(6>))total - («'(f))obs. (37) 

This uncertainty is plotted in Figure 2 (dotted line). For 
small survey radii, the variance in the reconstruction un- 
certainty is just the total variance of the convergence field, 
as the reconstruction uncertainty is dominated by missing 
structure beyond the survey boundary. As the survey radius 
approaches 9 = 0.1°, this effect begins to decrease, and be- 
yond 1 degree the survey is large enough to include all the 
relevant structure. 



[deg] 

Figure 3. The effect of a finite window on the lensing potential 
difference A</> for a Gaussian window of radius 9. The dotted line 
is the sampling variance including mode-mixing from the interior 
of a finite survey, while the solid line is variance due to missing 
modes exterior to the survey. Oscillations are real and due to the 
window function. The thick grey line is the shot-noise estimate 
for a survey with ri2 = 30 galaxies per sq. arcmin. 

only fully sampled for surveys above 10 degrees. On the 
other hand, when wishing to reconstruct cluster-scale mass 
concentrations, it is helpful to restrict a reconstruction to a 
cell-size of ~ 1 sq. deg. as this will cut out the large <j> fluc- 
tuations due to larger-scale structures (c.f. the good signal- 
to-noise obtained in this fashion for cluster reconstructions 
in Sections 5 and 6). 



3.2.4 Variance of the differential lensing potential 

We can also calculate the variance of the differential lensing 
potential field at the centre of a circular survey 



|A<M0)| 2 > 



d 2 £ 



r 4 c"|if<(o)i j 



1 - 



2Ji(ee) 



(£9) 



,(38) 



where Ji(x) is a Bessel function. The term in the square 
brackets subtracts off the mean field estimated over the sur- 
vey area. We plot this variance as a function of survey size, 
9, in Figure 3 (dotted line). We again assume a convergence 
power spectrum, C^ K , for a flat LCDM universe, and as- 
sumed the background galaxies are at z = 1. 

For small surveys the fluctuations expected in A<f> are 
greatly reduced, but on large scales the expected variation 
on A<f> becomes large; A(f> ~ 10~ 5 . This is due to the £~ 4 
weighting factor in equation (38) which makes the Acf) field 
sensitive to very large-scale structures. This is just a reflec- 
tion of the long-range nature of the 2-D lensing potential 
field, but means that the peak variance is on large angular 
scales. For smaller surveys these differential variations are 
suppressed as the survey becomes smaller than the struc- 
ture causing them. 

We conclude from this that the lensing potential field 
requires a large survey for a complete sampling of modes. 
For a LCDM model we find that the potential modes are 



3.2.5 Uncertainty in the differential lensing potential 

As well as the intrinsic variance of the differential 2-D lensing 
potential which we measure in a field, (|A</)(0)| 2 )ob s , we can 
also calculate the uncertainty in lensing potential due to 
missing modes from beyond the survey scale, (| A0(O)| 2 ) m is S : 



|A0(O)| 2 ) r 



|A<?i(0)| 2 ) t o t al - (|A^(0)| 2 ) obs 



(39) 



This is also plotted in Figure 3 (solid line). For small survey 
radii we again see that the uncertainty in the reconstruction 
is dominated by the missing shear structures beyond the 
survey boundary. In this case most of this missing structure 
is on larger scales. As we reach a survey radius of around 10 
deg. the survey begins to include this important large-scale 
shear structure, and the reconstruction uncertainty drops. 

Also plotted on Fig 3 is the shot-noise contribution to 
the uncertainty on a reconstruction from equation (20), as- 
suming 7i2 = 30 galaxies per sq. arcmin (thick solid line), 
typical for a ground-based survey. We see that shot noise 
is larger than the expected rms (j> fluctuations arising from 
large-scale structure for very small surveys (9 < 1 degree). 
On larger scales the shot noise is lower than the expected 
signal, allowing mapping of large-scale structure with good 
signal-to-noise on these scales. 

We note that the incompleteness contribution domi- 
nates over the shot-noise contribution for small surveys, 
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while the shot-noise contribution dominates for large sur- 
veys. The two contributions are around the same magnitude 
at around 9 = 15° when (A<j) 2 } w 2 x KT 11 . 

3.3 Wiener filtering lensing fields 

As we shall see, realistic galaxy ellipticities create a shot 
noise contribution to the gravitational potential reconstruc- 
tion which far exceeds the expected gravitational poten- 
tial amplitude from a cluster. This suggests one of two ap- 
proaches to reconstructing gravitational potentials in prac- 
tice: we can either examine the potential statistically from 
many objects of interest (e.g. stacking the signal from many 
groups or clusters); or we can filter the signal to overcome 
the large noise contribution. A valuable approach which we 
will use later involves Wiener filtering the gravitational po- 
tential (c.f. Hu & Keeton 2002). 

In order to apply this filtering, we constuct the vector P 
containing our gravitational potential measurements along 
a particular line of sight. We also calculate N, a matrix con- 
taining the noise covariance of the gravitational potential ra- 
dially along this line of sight; this can be measured directly 
from many $ reconstructions of a zero (j> field with appro- 
priate noise. Finally we require a matrix S, representing the 
expected covariance of the real gravitational potential signal 
along the line of sight. We use a multiple of the unit ma- 
trix for S, with amplitude chosen to equal the square of the 
expected gravitational potential amplitude, e.g. for clusters. 
Then we can apply a Wiener filtering 



R = S(N + S) _1 P 



(40) 



where R is our desired filtered gravitational potential. This 
filter uses our knowledge of the noise amplitude and covari- 
ance, together with the expected signal amplitude, to sig- 
nificantly reduce the impact of the noise. We will test the 
practicality of this approach in Sections 5 and 6. 

3.4 Photometric redshift errors 

In the above analysis, we assume that the distances to 
galaxies have been estimated from redshifts, measured ei- 
ther spectroscopically or photometrically. However, these 
redshifts include contributions from local velocities as well 
as the velocity of the overall Hubble flow. Thus the redshifts 
will cause a somewhat biased scatter in our distance mea- 
surements, as the local velocities are generated from gravita- 
tional instability due to the local Newtonian potential, and 
therefore correlate with our mass estimates. We must assess 
the level of error associated with this effect. 

We first consider the effect of a random distance er- 
ror, which is significant for photometric redshifts. We as- 
sume the source positions are perturbed by r — > r + e(r)r, 
where e(r) is a random field with zero mean and correlations 
(e(r)e(r')} = a 2 5o(r — r'). Expanding the lensing potential 
we find that the observed Newtonian potential field becomes 



<j>'{r) = 4>(r) + ed r <j){r). 



(41) 



This contributes to the first-order uncertainty in the New- 
tonian potential 



1 



r'$(r') 



(42) 



At large distances from the observer the terms in the square 
brackets vanish, and the leading contribution to the uncer- 
tainty in distance comes from the gradient of the Newtonian 
potential. Since <r e is, for appropriate redshift surveys, small 
in comparison with the redshift depth probed, and since 9 r $ 
will be small in a smoothed survey, this effect should not be 
dominant in recovering the gravitational potential. To con- 
firm this, we will examine the effect of redshift errors on our 
simulations in Section 6.1. 



3.5 Redshift-space distortions 

We must now consider the effect of velocity distortions, 
which may be significant for the more accurate spectroscopic 
redshifts. In linear theory, the velocity field is related to per- 
turbations in the mass-density field by 



v = -Hf(Q m ) v\7-' 2 s 



(43) 



where f(0. m ) = dlnS/dlna w Q,^ 6 is the growth index of 
density perturbations (e.g. Peebles 1980). The position of 
galaxies are then shifted into redshift space by 



r — > s — r + u(r)r 



(44) 



where s is the redshifted position and u(r) = f.v(r) is 
the radial component of the velocity field. The distorted 
redshift-space lensing potential is then 



(45) 



In this case the systematic distortion of the Newtonian po- 
tential is second-order: 



A<E> = 



/ 2oA 



H 1 [k$$ + 2$^$ + (<9 r $) 2 ] 



(46) 



The magnitude of this effect on the reconstructed density 
field is A<5 w fS 2 , and so only contributes to second order. 



3.6 Multiple scatterings 

A further concern is the fact that a fraction of light rays will 
be multiply scattered as they travel from source to observer. 
How will this affect our assumption that the shear field can 
be derived from a lensing potential? 

The scattering of light rays can be written as 



91 = D ij {e)59 j 



(47) 



where Dij is the lens distortion matrix, defined for a single 
scattering by 



Da = 8« + didj<t> = (1 - + la- 



(48) 



The distortion matrix for n multiple scatterings is then just 
the product of distortion matrices, 



D n = D W fci D (2) fc 2 _ _ _ D („) 

IJ % fci kn] 



(49) 



where is the effect of scattering off the i th structure 
along the light path. 

For the case of double scattering we can then define an 
effective convergence, shear and a rotation, 
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toff = m + K-2 — 2kik 2 — Tr 7j ^72 kj 

lif = 7l ij + 72 ij - ^172 ij - K271 ij 

+7i fc (i72j)fc - [Tr7! fc ^ fe ]5^ 

Wij = 7i fc [i72fc 3 ]- (50) 

Thus, even for a double scattering, the presence of a rota- 
tional component to the distortion matrix shows that the 
distortion can no longer be strictly constructed from a po- 
tential. However, in the case of weak lensing at the 10% 
level, we see that the rotational components will on average 
be around 0.1%. However we can expect this to break down 
in the strong lensing regime, very close to massive cluster 
centres. Therefore it is worth considering how common pro- 
jections may be in lensing. 

If we assume that the cluster spatial distribution is ran- 
dom, the probability of finding one or more clusters at ran- 
dom along a given line of sight of volume V = 4ty/R 3 /S 
is 

P(N > 1) = 1 - e- x « A « N c f, (51) 

where N c is the number of clusters in the whole sky volume 
and / is the angular fraction of the sky covered by our line 
of sight. The probability of seeing two or more clusters along 
a given line of sight is 

P(N > 2) = \\ 2 (52) 

The statistic we require is the conditional probablity of find- 
ing a second cluster when we have found one, 

This yields 

where 6 C \ is the angular size of a typical cluster. This calcula- 
tion is highly approximate, but allows us to see that finding 
clusters behind clusters can occur with non-negligible prob- 
ability. This will cause no difficulties for our method, except 
in the strong lensing regime very close to the centre of the 
clusters. Indeed, our method is a useful means of measuring 
several mass concentrations along the same line of sight. 



4 SIMULATING 3-D LENSING 

In order to investigate the practical application of the re- 
construction formalism, we have conducted a series of sim- 
ulations representing a realistic space volume for a lensing 
survey, including galaxies with an appropriate spatial dis- 
tribution and intrinsic ellipticity and gravitational lenses of 
appropriate size and mass. The shear and magnification for 
objects behind the lenses can be calculated (retaining the 
full information regarding 3-D variation of these quantities), 
and the objects' shapes can be altered accordingly. With this 
flexible simulation package in place we can attempt to recon- 
struct the gravitational potential causing the lensing, using 
knowledge of only the galaxy ellipticities, their redshifts, and 
equations (8) and (9). Here we describe in detail the form 
of these simulations. 



Immediately we are faced with a question as to which 
of the fields described so far (7, k, <f>, 8) we should use for 
3-dimensional analysis. In reality, the most appropriate field 
to use depends upon the application intended. For detection 
and measurement of mass concentrations along the line of 
sight (Section 7), the fields k and <j> are the most useful, as 
they have the best signal-to-noise (~ 7a in each case in [9 
arcmin 2 , Az — 0.05] pixels, for field radius 0.5°; see Section 
6) and are transversely local representations of the mass 
present. For a direct mapping of the gravitational field, $ 
is most appropriate for mapping particular concentrations 
such as groups and clusters, as it has much less noise than 
the 8 field on < 1° scales (c.f. equation (25)). Indeed, as 
the noise in the 8 field grows quadratically with survey size 
(equation (27)), as does the noise in the $ field, it may be 
that <I> is most useful at large survey areas as well, depending 
on the required application. In this paper we concentrate on 
reconstructions of cluster size mass concentrations, and will 
therefore make use of the <j> and $ fields. 



4.1 Constructing the Shear Field 

We construct a 3-dimensional grid, typically with a total 
of 100 3 points. This represents the redshift cone in which 
we will attempt to reconstruct the gravitational potential, 
i.e. the x and y directions represent an angular range on 
the sky, while the 2 direction represents a radial distance. 
Typically we will use this to model angular scales in x and 
y of 1° while probing in the z direction down to an effective 
redshift of 1. We will therefore quote coordinates in r = 
(x,y,z) ranging from to 1. This does not imply that the 
z coordinate represents redshift, however; throughout this 
paper, z represents a comoving distance measure. 

This choice of coordinates significantly simplifies our 
analysis: in a flat universe, the comoving transverse separa- 
tion of unperturbed light ray paths converging at an observer 
is proportional to the comoving radial distance along the 
light path (see e.g. Bartelmann & Schneider 2000, Section 
6.1). Thus unperturbed rays will simply move along our 2 
coordinate with fixed (x,y). Also, in order to lay down e.g. 
a constant 3-dimensional comoving number density of ob- 
jects, we simply allocate a constant number density in our 
coordinates; no correction is necessary for a varying physical 
number density. Finally, the unitless transverse derivatives 
required in equation (2) are simply d x and d y in our coordi- 
nate system. 

We fix in this grid the positions of lenses which we will 
wish to recover. We assign to each lens a mass and three 
perpendicular scale lengths for the size of its gravitational 
potential. In the following section we will be concerned with 
typical galaxy cluster lenses, in which case we assign masses 
in the range 0.5-5X 1O 14 M and radii of 0.5-2Mpc. We will 
position these clusters between z — 0.25 and 2 = 0.5 while 
probing lensed galaxies down to 2 = 1, mimicking lensing 
studies of massive clusters (e.g. Tyson et al 1990, Kaiser & 
Squires 1993, Bonnet et al 1994, Squires et al 1996, Hoekstra 
et al 1998, Luppino & Kaiser 1997, Gray et al 2002). 

Given these masses and scale lengths, we can calcu- 
late the gravitational potential over the entire 3-dimensional 
grid. We have used Navarro, Frenk & White (1996) profiles 
for the density, 
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Figure 4. Example of gravitational potential. Here we simu- 
late the gravitational potential for two NFW clusters at r = 
(0.5,0.5,0.25) and (0.5,0.7,0.4); the upper panel is a 3-D rep- 
resentation of (y, z,<&) resulting from an x slice at x = 0.5°; the 
lower panel displays the $ values for (y, z) in greyscale. 

p(r) =q/\r-r c \(l + \r-r c \) 2 , (55) 

where q is a measure of the mass, and r c are the coordi- 
nates of the centre of the mass profile. Using Gauss's law 
and a further integration and renormalisation, we find the 
resulting gravitational potential for a mass concentration 

*(r) = -^Iog(l + J^), (56) 

with qi as a measure of mass and r s the NFW scale pa- 
rameter. The total gravitational potential in the simulated 
volume of space is taken to be the sum of the individual lens 
contributions. 

Figure 4 shows an example of our constructed gravi- 
tational potential, for two NFW profile clusters, placed at 
z — 0.25 and z — 0.4; the view is of a 1 grid-unit slice at 



x — 0.5°. The apparent distortion of the clusters is due to 
the scales chosen; we are examining a large distance scale 
radially (~ 3000Mpc), with much smaller distance scales 
transversely (~ lOMpc). We will wish to reconstruct this 
gravitational potential for all (x, y, z) from shear informa- 
tion only. 

The masses of the clusters in this example were chosen 
to mimic the shear properties of clusters in Gray et al (2002), 
with mass m = 8 x 1O 13 M0 within a radius of 2 arcmin. 
Note that in order to calculate <f> and $ we convert the x, y 
coordinates from degrees to radians, leading to the same (f> 
and <J> normalisations as in Section 3. 

From the gravitational potential shown in Figure 4, we 
calculate the lensing potential (j> given by equation (7) . This 
is a necessary step towards calculating the shear, which is 
all we will use for potential reconstruction. For each point in 
(x, y) we set 4>(x, y, 0) = and then use the discrete version 
of equation (7), 

z 

<t>(x,y,z) = -2^2 ^-^-$(x,y,w)Aw, (57) 

w = l 

where we must first use this equation to calculate 
4>(x, y, 0.01) then <f>{x, y, 0.02), etc. In this fashion, we can 
calculate the lensing potential for all points on our 3- 
dimensional grid. 

Figure 5 shows the lensing potential calculated as above 
for the example introduced in Figure 4. Note that the axes 
on this 3-D plot are (y, z, <f>); we are observing how the lens- 
ing signal grows with depth. The figure shows the generic 
behaviour for all lensing; the lensing potential due to distor- 
tion from a massive object becomes stronger with increasing 
depth, but asymptotes to a finite value at large z. 

We are now in a position to calculate the 3-dimensional 
gravitational shear field arising from the gravitational po- 
tential. We calculate the shear components 71 and 72 from 
our 3-dimensional 4> field by first approximating 

9 xx </>(x) ~ [4>(x + Ax,y,z) + 4>(x - Ax,y,z) 

-2<j>(x,y,z)]/(Ax) 2 . (58) 

An entirely similar approximation is made for d yy (f)(x), while 
we approximate 

d xy 4>(x.) ~ ^[()>(x + Ax,y + Ax,z) + 4>(x — Ax,y — Ax,z) 

-4>{x - Ax, y + Ax, z) (59) 
-<f>(x + Ax,y~ Ax,z)]/(Axf. (60) 

Then, following from equation (2), we can write the shear 
components as 

7i (x) = -i(cW(x) - d yy <t>(x)) (61) 
72 (x) = -d xy 4>(x) (62) 

From these equations we can calculate shear values for all 
points on our 3-dimensional grid. An example of a calcu- 
lated shear field is shown in Figure 6, corresponding to the 
gravitational potential of Figure 4; this is an (x, y) slice of 
the 3-D shear field at z — 0.75. Note the clear signatures 
due to the two clusters. A slice further back in z would have 
a similar x,y pattern but a larger shear magnitude; a slice 
in 2 in front of the clusters would have no shear signal. 

We have therefore calculated the gravitational shear 
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Figure 5. Lensing potential calculated for the gravitational po- 
tential of Figure 4. Note the increase in distortion expected with 
growing z. In the upper panel, the axes represent the (y, z,(p) co- 
ordinates; in the lower panel the <j> values for (y, z) are shown in 
greyscalc. 



which would exist for an object at any grid point; now we 
must place objects at some of these grid positions, with an 
appropriate distribution in space. We normalise the total 
number of objects to those expected in surveys down to the 
required depth; for a ground-based survey probing to a lim- 
iting redshift z r = 1, we expect a number density n ~ 30 
usable galaxies per square arcmin (e.g. Bacon et al 2002). 
Alternatively, we can project the number densities expected 
for deep, space-based surveys down to z r — 1 (n ~ 100; see 
Massey et al 2002). 

The selected total number density is used to find the 
probability of an object existing at a given grid point in the 
following fashion. We treat each (x, y) sheet as a slice of a 
comoving pyramid (not a box), in order to take into account 
the fact that the survey has an angular extent. In some of the 



1.0 
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Figure 6. Slice of the 3-D shear field corresponding to the grav- 
itational potential of Figure 4. This is an (x, y) slice at z = 0.75. 
The largest 7 value plotted is 0.27. 

simulations we will assume that the 3-dimensional number 
density is constant, while in others we will adopt dN/dz oc 
z 2 exp(— az 2 ); the distribution will be stated in each case. 

In the case of varying N(z), the number of galaxies N 
expected per grid point can be found using 

KTI . dN Az ,„„. 

N ^ = ^w (63) 

where 9 is the transverse survey length (usually 1° in our 
simulations), Az is the radial grid increment, and N to t = 
nO, where O is the solid angle extent of the survey. For each 
z, we then generate N(z)9 2 random coordinates in x and 
y, and increment the number of objects at the nearest grid 
point for each coordinate pair. 

For each object, a random Gaussian-distributed shear 
value is chosen to simulate the effect of intrinsic ellipticity, 
and is added in quadrature to the gravitational shear value. 
For typical ground-based surveys the scatter in shear esti- 
mators due to the intrinsic ellipticity can be well modelled 
by Gaussians in 71 and 72 with standard deviation 0.3 in 
each component (e.g. Bacon et al 2002). For space-based 
survey simulations, we use the smaller standard deviation 
of 0.2 in each shear component (see e.g. Rhodes et al 2001). 

We now have a set of galaxies sampling the 3-D shear 
field at a finite set of points in space, with an additional 
shot noise contribution from their intrinsic shapes. We must 
now attempt to estimate the underlying purely gravitational 
3-D shear field, and from that to calculate the 3-D lensing 
potential 4> from equation (9), and the corresponding 3-D 
gravitational potential $ from equation (8). 

4.2 Reconstructing the Potential 

We can overcome the shot noise from galaxy ellipticities by 
a combination of binning many galaxies' shears in a cell, 
and smoothing the shear field. In the specific examples given 
below, we find it convenient to initially rebin our galaxy grid 
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to pixels with 3 arcmin diameter in the transverse direction 
and Az = 0.05 in the radial direction; further smoothing can 
be carried out as necessary later. We calculate the averaged 
shear field for the rebinned grid; this averaged shear will 
then be the basis of our potential reconstruction. 

We can conveniently find the <f> field corresponding to 
this 7 field using Fourier tranforms of these two fields. Fol- 
lowing Kaiser and Squires (1993), we find from equation (9) 
an optimal estimate of the lensing potential, 



0(k) 



2{kl -^)7i(k)+4fefc a72 (k) 



(64) 



where </>(k) etc are Fourier transformed quantities. This 
equation is valid if the transverse differential operators of 
equation (9) are dimensionless; this is satisfied by our model, 
where the x and y directions represent angles rather than 
physical distances. 

In this fashion we calculate 0(k) using Fast Fourier 
Transforms, and hence find 4>- Since noisy low-fc modes can 
introduce rather large variability in with z (see Section 
2), we further add a constant at the edge of our (real-space) 
field to ensure that the mean of (j> around the edge of the 
field is zero (this is equivalent to equation (12) above; see 
Section 3.1.2 and 3.2 for a discussion of the level of error this 
causes). We can then smooth the <f> field in the z direction 
to reduce the noise amplitude, convolving the <f> field with a 
radial Gaussian. We can choose the width of this kernel in 
each circumstance as a compromise between reducing noise 
and retaining spatial resolution. 

We then calculate the gravitational potential $ using 
equation (8). As in equations (58) and (60), we use local 
differences in <f> in the z direction to approximate the deriva- 
tives. If necessary, we can smooth or filter $ itself to further 
reduce noise levels. 



5 SIMULATION RESULTS 

We have described above our means of simulating mass dis- 
tributions and corresponding shears, followed by our method 
for reconstructing the lensing potential and gravitational po- 
tential given a finite number of noisy estimators of the shear. 
Here we describe our findings for this reconstruction. 




5.0x10 8 1.0x10 7 1.5x10 7 2.0x10 7 2.5x10 7 3.0x10 7 3.5x10 7 




5.0x10 J 1.0x10 u 1.5x10 u 2.0x10 



5.1 Perfect reconstruction 

First we examine the accuracy of our method when the shear 
field is perfectly known everywhere. Figure 7 shows the re- 
construction of the lensing potential when we simply input 
into our reconstruction the full shear field shown in Figure 
6; we do not smooth the shear field in this case. We regain 
the lensing potential successfully; this is quantified in Figure 
7c, where we show the difference between the reconstructed 
and original <j) fields, for a slice through the reconstruction. 
This error introduced by our numerical implementation is 
< 3% of the lensing signal everywhere behind the clusters 
within 0.15 degrees of the cluster centres. The one exception 
is the core pixel line behind each cluster, where the error is 
6.6% and 9.7% of the lensing potential for the left and right 
cluster respectively. This is expected, due to the cusp within 
this pixel, which the binned shear cannot accurately follow. 
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Figure 7. Top and middle panel: reconstructed lensing potential 
using the full shear field of Figure 6, in a 100 3 grid. This is a 
(y, z) plane at x = 0.5°. Bottom panel: difference between input 
and recovered lensing potential fields. 
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Figure 8. Top panel: reconstructed gravitational potential using 
the full shear field of Figure 6, in a 100 3 grid. This is a (y, z) 
plane at x = 0.5°. Bottom panel: difference between input and 
recovered gravitational potential fields. 



We see that our procedure for reconstructing the lensing 
potential is very successful in the absence of noise. 

We proceed to use our reconstructed lensing potential to 
find the 3-D gravitational potential, still using the full shear 
field (i.e. still without sampling by a finite number of objects, 
or adding shot noise due to galaxy ellipticities). Our result is 
shown in Figure 8, which displays the reconstructed $ and 
the difference between input and reconstructed fields. Again 
it can be seen that, with perfect knowledge of the shear 
field, a good reconstruction is achieved with our method. 
For the cluster at z = 0.25, the error is < 0.4% of the signal 
within a radius of 0.2 degrees of the cluster, except for the 
core pixel where the error is 11.5%. For the cluster at z — 
0.5, the error is < 5% of the signal within a radius of 0.1 
degrees of the cluster, except for the core pixel where the 
error is 17.6%. This is again due to the cusp at the cluster 



centre, which cannot be followed well by our averaged shear 
field. Nevertheless, it is clear that our method is successful in 
reconstructing cluster gravitational potentials in the absence 
of noise. 



5.2 Reconstruction with noise 

Having demonstrated that the inversions of the shear field 
to obtain lensing and gravitational potentials are viable in 
the absence of noise, we now wish to add the two primary 
sources of noise present in lensing experiments: Poisson noise 
due to only sampling the field at a finite set of galaxy po- 
sitions, and additional noise due to galaxies' non-zero ellip- 
ticities. We use appropriate number densities for plausible 
ground-based experiments (30 per sq arcmin) and space- 
based experiments (100 per sq arcmin), and use equation 
(63) to place an appropriate number of objects at each red- 
shift slice. We continue to use the gravitational potential of 
Figure 4. 

We incorporate the effect of the intrinsic ellipticities of 
the galaxies as described in Section 4 (i.e. adding a Gaussian- 
distributed random shear value to the gravitational shear, 
with a standard deviation of 0.2 per shear component for 
space-based applications and 0.3 for ground-based applica- 
tions). 

Given this noisy shear field, we carry out our reconstruc- 
tion as described in Section 4, with a Gaussian smoothing 
of the <f> field in the z direction with a la width of 0.1 in 
redshift (not applied for calculating Wiener-filtered <E>). 

We will now discuss the reconstructions obtained for 
ground-based and space-based data. In the discussion below, 
we will often describe measurement significance in terms of 
f(x) = J(x)/cr(x), with the amplitude of the lensing 
potential or gravitational potential at a particular point x, 
and cr(x) the noise level in the vicinity of this point. 



5.2.1 Lensing Potential for Space-based Experiment 

The resulting lensing potential reconstruction for our space- 
based experiment is shown in the top panel of Figure 9. 
(In this and later figures, we plot the results for our 20 3 
grid, but display a resampled grid calculated by padding the 
Fourier transform of the grid with high-fc modes set to zero.) 
We find that we obtain a reasonable reconstruction of the 
lensing potential, with v ~ 6.9 per pixel in the background 
(z > 0.75) for the nearer cluster, with v ~ 4.2 for the z = 
0.4 cluster, within 0.1 deg radius of the cluster centres. We 
could increase this signal by rebinning or smoothing, at the 
cost of reducing spatial resolution. We could also find an 
overall signal-to-noise for each cluster by finding a means of 
averaging all of the lensing signal arising behind a cluster; 
we will discuss this in Section 7. Note the large noise peaks 
in the foreground (z < 0.2) of the reconstruction, due to the 
small number of objects available in this volume. 

The lower panel of Figure 9 shows the difference be- 
tween the input and recovered lensing potential. Pleasingly, 
we observe no evidence of residuals associated with miscon- 
struction of the lensing potential. The noise levels are as 
expected from equation (20), as discussed in Section 6 be- 
low. 
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Figure 9. Top panel: reconstructed lensing potential A</> using 
finite number of galaxies with realistic ellipticities; n = 100 per 
sq arcmin, 7 rm s = 0.2 per shear component, as expected for a 
notional space-based survey. Crosses show the positions of the 
cluster centres. Bottom panel: difference between input and re- 
covered lensing potential fields. 

5.2.2 Gravitational Potential for Space-based Experiment 

On the other hand, we find that the full 3-D reconstruction 
of the gravitational potential for typical cluster masses is dif- 
ficult even from space unless we include filtering (see Section 
6.1). For example, a measurement of the gravitational po- 
tential for a set of 25 clusters with mass m = 8 x 10 13 within 
2 arcmin added together at z = 0.25, having smoothed (j> in 
the radial direction with a top hat with Az — 0.25, produces 
a gravitational potential amplitude for the central pixel of 
the cluster which is only 3.1 times larger than the rms noise. 
However, the noise amplitude for larger z grows rapidly (c.f. 
Section 6.1), making 3-D measurement of a resolved cluster 
potential difficult without filtering. 

In order to improve on this, we can use the Wiener 
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Figure 10. Reconstructed gravitational potential for our space- 
based experiment, using galaxy properties as in Figure 9, af- 
ter Wiener filtering. Note the detections of clusters at (y, z) = 
(0.5,0.2) and (0.7,0.4) with u = 4.2 and 2.1 respectively (crosses 
show trough minima). 
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Figure 11. Reconstructed gravitational potential for a notional 
ground-based survey, after Wiener filtering; n = 30 per sq arcmin, 
7rms = 0.3 per shear component. The clusters are measured with 
v = 1.9 and v = 1.4. 



filtering described in Section 3. We create a vector P con- 
taining our measured gravitational potential along each line 
in the z direction, and measure the covariance of the <3? noise 
along this line of sight from 100 zero-0 reconstructions of $, 
recording this in a matrix N. We set the signal covariance 
matrix S = (3 x 10 - ) 1 where 1 is the unit matrix, in order 
to select for a signal expected for a small cluster mass (cf 
Figure 8) . We then apply equation (40) to our gravitational 
potential vector for each line of sight. 
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Figure 10 shows the gravitational potential measure- 
ments for our notional space-based survey after Wiener fil- 
tering (this is again a reconstruction for the single field 
shown in Figure 8; after Wiener filtering we no longer need 
to stack many fields to obtain a signal). We measure the 
z = 0.25 cluster with v — 4.2 at the peak (trough) pixel of 
its gravitational potential well, with v — 2.1 at the gravita- 
tional potential trough of the z = 0.4 cluster. We also find 
a substantial noise peak in the foreground at z = 0.1. While 
the recovery of the cluster gravitational potentials therefore 
constitutes a challenging measurement, we are indeed able 
to reconstruct useful information in the gravitational poten- 
tial field itself. (The detection significance of these clusters 
is much higher than the measurement v at a given point in 
the cluster; see Section 7 for an approach to the detection 
significance.) Note the reduced absolute amplitude of the 
gravitational potential in Figure 8; this is due to the scaling 
in equation (40). 

5.2.3 Gravitational Potential for Ground-based 
Experiment 

From the ground, we find that we can again reconstruct 
the lensing potential, with v ~ 2.3 for 0.75 < z < 1.0 
given a cell size of 0.05 in z. As before, we could improve 
the signal-to-noise by reducing our spatial resolution. How- 
ever, reconstruction of the gravitational potential itself (with 
Wiener filtering) is more challenging than from space: Fig- 
ure 11 demonstrates that we recover the z — 0.25 cluster 
with v — 1.9 at its centre pixel, but have a less significant 
measurement of the second cluster amplitude with v = 1.4 
at its centre. The second cluster has the expected position in 
(x, y), but is offset to z = 0.3; this degree of offset, Az = 0.1, 
is found to be typical for ground-based attempts at measur- 
ing the 3D gravitational potential, due to the high $ noise 
level making Wiener filtering somewhat inaccurate in the z 
direction. As with our space-based experiment, we also find 
a substantial noise peak in the foreground at z = 0.1. 



6 PROSPECTS FOR MAPPING 

The results above are encouraging for the mapping of the 
3-D lensing fields. This includes the gravitational potential; 
we can use Wiener filtering to detect individual mass con- 
centrations, or can stack the noisy potential field from many 
clusters in order to obtain information on the typical grav- 
itational profile of mass concentrations. Here we will exam- 
ine the prospects for mapping with the various 3-D fields 
we have discussed so far, paying close attention to the noise 
contributions to each field. 

6.1 Noise amplitudes 

In order to demonstrate the feasibility of our reconstruc- 
tions, we examine the noise from convergence, lensing po- 
tential and gravitational potential maps measured in our 
simulations as a function of z. Figure 12 shows the rms noise 
amplitude for function of redshift; Figure 13 shows 

the corresponding uncertainty in <j>, while Figure 14 shows 
the uncertainty in $. In each case, we show the measured 




z 



Figure 12. Noise variance for 1 pixel (0.05 deg, 0.05 in z) as a 
function of z for the 3-D convergence field k. The thin solid line 
represents the measured noise level for our fiducial space-based 
experiment, while the dashed line represents measured noise from 
our ground-based experiment. The dotted line shows the expected 
signal from a cluster at z = 0.2. The thick solid line shows the 
noise expected from theory for the ground-based survey, from 
equation (15); we see that there is close agreement between sim- 
ulations and theory. 
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Figure 13. Noise variance for 1 pixel (0.05 deg, 0.05 in z) as a 
function of z for the lensing potential in a 1 deg 2 square survey. 
The solid line represents the measured noise level for our fiducial 
space-based experiment, while the dashed line represents mea- 
sured noise from our ground-based experiment. The dash-dotted 
line shows the effect of redistributing the galaxies according to 
equation (63) for our space-based experiment. The dotted line 
shows the expected signal from a cluster at z = 0.2. The thick 
solid line shows the noise expected from theory for the ground- 
based survey, from equation (20); we again find good agreement 
between simulations and theory. 
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Figure 14. Noise variance for 1 pixel (0.05 deg, 0.05 in z) as a 
function of z for the gravitational potential. The solid line rep- 
resents measured noise level for our fiducial space-based experi- 
ment, while the dashed line represents measured noise from our 
ground-based experiment. Top: standard gravitational potential 
reconstruction. The dotted line shows the expected signal from a 
cluster at z = 0.2. The thick solid line shows the noise expected 
from theory for the space- based survey, from equation (22); we 
again find good agreement between theory and simulation. Bot- 
tom: Wiener filtered reconstruction. The spikes show the mea- 
sured Wiener-filtered amplitude of a 1O 14 M0 cluster at z = 0.2 
and z = 0.4 for ground-based (dotted) and space-based (solid) 
noise levels. 



noise amplitude for space-based and ground-based experi- 
ments, measured as a function of z, along with the expected 
potential amplitude from a m = 10 14 M Q cluster. Here we 
have used the averaging scales described above, i.e. we have 
calculated the noisy k and 4> fields on a grid with total width 
1 degree, with grid spacing of 0.05 degree and 0.05 c/Ho in 
z. It is simple to derive noise amplitudes for other survey 
configurations from equations (19) and (22). 



We see from from the theoretical curves on Figure 12 
that the noise amplitude measured for the convergence field 
is in good agreement with that calculated in Section 2. Sim- 
ilarly, we see that the measured noise amplitude for cf> in 
Figure 13 agrees well with the predicted theoretical curve. 
Lastly, we see in Figure 14 that the uncertainty in gravita- 
tional potential is in agreement with our theoretical model. 

Note that, for the convergence and lensing potential, 
the noise amplitude is such that v ~ 6 measurement of these 
fields from a typical z = 0.2 cluster is possible in (0.05 de- 
gree, 0.05 in redshift) bins in a space-based experiment, if we 
examine the background potential at z > 0.5 in the centre 
of the cluster. For a ground-based experiment, the signal-to- 
noise in these bins is ~ 3. Thus as we saw above, it is possible 
to map the properties of clusters in terms of the lensing po- 
tential from space or ground. Note the effect in Figure 13 
of redistributing the galaxy distribution according to equa- 
tion (63); the amplitudes of the noise are comparable, with 
a slight decrease in noise at low redshift and a correspond- 
ing increase at high redshift for the galaxy distribution of 
equation (63). 

On the other hand, the gravitational potential without 
Wiener filtering is only measured at v — 0.38 in equiva- 
lent bins around z = 0.2 for space-based experiments, and 
has a signal-to-noise of only ~ 0.12 for ground-based mea- 
surements. Since the noise reduces as V^V where N is the 
number of stacked fields, we would have to stack ~ 40 fields 
in order to resolve the gravitational potential at the v = 2.5 
level, for a space-based experiment without Wiener filtering. 
Increasing the size of bins is not an option in this case, as 
we will lose spatial resolution for examining the profile of 
the cluster. 

However, Wiener filtering allows us to recover cluster 
masses successfully, as shown in the bottom panel of Figure 
14. Here we see that, for a cluster mass of 10 14 M Q at z = 

0. 25, we obtain v — 2.1 measurements from the ground and 
v = 3.3 measurements from space. Note the differing noise 
and signal levels expected after Wiener filtering for ground 
and space-based experiments; this is due to the differing 
weightings in equation (40) given different input noise levels. 

6.2 Dependence on transverse and radial pixel 
sizes 

The effect of increasing the angular pixel scale in a survey 
is of interest, as this might be thought to increase signal- 
to-noise. However, Figure 15 shows how an increase in pixel 
size for a fixed survey size reduces the lensing potential noise 
variance. We note that the effect is small until a pixel size 
~ 0.25Q 3urvey is used; at this point, the pixel size will usu- 
ally be far too large to be of use to us, as we will typically 
be interested in spatially resolving objects with the lensing 
potential. The cause of the decrease in noise as pixel size in- 
creases is partially the In r force law involved in the analysis 
in Section 3.1.2, and partially the increased bin size leading 
to averaging out of the small-scale noise. 

Figure 16 shows the effect of increasing the radial bin 
size for the gravitational potential reconstruction. We see a 
reduction of the noise level in agreement with equation (22) , 

1. e. cr($) oc (Az)- 5/2 . Thus if we are unconcerned with ra- 
dial resolution, we can increase our signal-to-noise for grav- 
itational potential by increasing Az; however, we will often 
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Figure 15. Noise variance for the lensing potential at z=l as 
a function of angular pixel size, for our space-based experiment, 
with a 1 deg 2 square survey. 
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Figure 17. Top panel: reconstructed lensing potential using the 
full shear field of Figure 6, while including redshift uncertainty 
Az = 0.1. This is the usual (y, z) plane at x = 0.5. Bottom panel: 
difference between input and recovered lensing potential fields. 



Figure 16. Noise variance for the gravitational potential as a 
function of pixel size in the radial direction, for our space-based 
experiment (at z=0.4, 1 deg 2 survey, angular pixel diameter 3'). 
The dotted line is the expected dependence from equation (22). 



be attempting to locate a mass concentration in the radial 
direction, so this procedure should be used with caution. 



6.3 Redshift Errors 

In the above simulations, we have demonstrated that the 
Poisson noise of sampling from a finite set of galaxies, to- 
gether with the noise due to galaxy ellipticities, represent 
serious sources of uncertainty for our reconstruction, which 
must be overcome by averaging the signal within sufficiently 
large pixels or filtering the signal. However, there remains 



the further source of error due to redshift measurement un- 
certainty, which we examine here with our simulations (see 
Sections 3.4 and 3.5 for analytical discussion). 

We can examine redshift errors with our simulations by 
allowing an uncertainty in the redshift of each galaxy, in 
the following fashion. For each galaxy, a z value is drawn 
uniformly within the z slice in question; an uncertainty is 
introduced in this z coordinate, by drawing from a Gaussian 
random variable with la width Az = 0.1 (pessimistically, for 
photometric redshifts; cf Brown et al 2002 with Az — 0.05 in 
< z < 0.8). If the new z coordinate (resulting from adding 
this random offset to the original z position) is moved to a 
new shell, the galaxy (with its shear calculated for the slice 
which it intially belonged to) is moved to the neighbouring 
slice. 

Figures 17 and 18 show the result of this process when 
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Figure 18. Reconstructed gravitational potential using the full 
shear field of Figure 6, including redshift uncertainty Az = 0.1. 

the shear is fully known everywhere; note that even with 
this large redshift uncertainty, the change in the lensing po- 
tential reconstruction is small (~ 10% at redshifts near the 
cluster redshift, and less elsewhere). Note that this smearing 
effect leads to slightly higher lensing potentials in front of 
the cluster, and slightly slower rise in the potential behind 
the cluster. Thus, for lensing potential reconstructions, this 
effect will not dominate the noise. However, Figure 18 shows 
that the result of such an uncertainty will be a smearing 
of the cluster gravitational potential with smearing width 
~ Az; understandably, we cannot reconstruct the potential 
with a resolution greater than our redshift resolution. 



7 3-D INFORMATION FROM THE LENSING 
POTENTIAL 

It will be noted from our simulations that recovery of the 
gravitational potential is more difficult than an adequate re- 
construction of the lensing potential. It is easy to see why 
this should be so; the gravitational potential requires dou- 
ble differentiation of the already noisy lensing potential field. 
Therefore we are interested in both the lensing and gravita- 
tional potentials; while the lensing potential is more easily 
accessible, and itself contains useful topological information 
about the mass field, the gravitational potential is the more 
fundamental quantity. 

A particular example of what is achievable by studying 
the 4> field is the characterisation of the 3-D matter distri- 
bution on cluster scales. Treating a cluster as a mass delta 
function in the z direction (c.f. Hu & Keeton 2002), it is 
clear from equation (7) that the expected lensing potential 
in a flat universe due to a cluster at radial position r 3 will 
be 

= { IpdV^f r > I] ^ 
where pdV is the mass content of the source pixel. 
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Figure 19. 1 and 2<x \ 2 fit constraints on z and m for 1 cluster 
along a line of sight, m = 8 X 1O 13 M0, with noise appropriate for 
our space-based experiment. Here we are fitting the simulated 3-D 
<j> field with the <j> field from a 2-D NFW profile with parameters 
(z,m). 
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Figure 20. Best fit model for (f> field along line of sight through 
centre of a cluster with m = 8 X 10 14 Mq, with noise appropriate 
for our space-based experiment. 



We can use a superposition of such cluster contributions 
to fit a given <f> field, and thus discover significant clusters 
behind clusters, for example (c.f. Hu & Keeton 2002). We 
can also directly measure constraints on mass and position 
of clusters without using the redshift of the cluster members 
themselves (e.g. Wittman et al 2001, 2002). 
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Figure 21. 1 and 2a x 2 fit constraints on zi and Z2 for 2 clus- 
ters along line of sight, with noise appropriate for our space-based 
experiment, having marginalised over the masses of the two clus- 
ters. 

7.1 Single cluster \ fitting 

In order to demonstrate these applications, we first simu- 
late a cluster at z = 0.25 with mass 8 x 10 13 Mq within a 2 
arcmin radius using the simulation recipe described in Sec- 
tion 3, including realistic galaxy distribution and ellipticity 
(using our space-based parameters, n — 100, a 7 = 0.2 per 
component). After measuring the resulting <f> field as in Sec- 
tion 4, we applied a \ 2 fitting procedure for the mass m and 
position z of the cluster. This was achieved by setting up a 
2-D NFW profile at radial position z, normalised to mass 
m; the expected 3-D 0-field for this profile was calculated 
according to equation (65). \ 2 f° r the data with respect to 
this profile was calculated for an array of m and z values, 
with step size 2.5 x 10 Mq in m and 0.005 in z. Note that 
the x, y position and radius of our NFW test profile were 
fixed to the actual position and radius of the cluster in this 
experiment; in a practical scenario one would apply a \ 2 fit 
for these parameters as well, or infer them from the galaxy 
positions of cluster members. 

Figure 19 shows the resulting constraints on mass and 
radial position; we find a highly significant detection of the 
cluster at the Ax 2 = 186 level, thus we can certainly use 
this 3-D approach to detect at least one cluster along the 
line of sight. We obtain accurate measurements of the mass 
(m = (8.07 ±0.83) x 10 13 M Q within 2' radius) and position 
(z — 0.24 ±0.016), which makes this approach promising for 
examining clusters in 3-D. Figure 20 demonstrates a best-fit 
(f> field for this form of simulation, with mass multiplied by 
10 for illustrative purposes. 

7.2 \ 2 fitting for two clusters along the line of 
sight 

We can examine the possibility of detecting clusters behind 
clusters by simulating an 8 x 10 13 M© cluster (2 arcmin ra- 
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Figure 22. 1 and 2cr \ 2 fit constraints on mi and m,2 for 2 
clusters along line of sight, with noise appropriate for our space- 
based experiment, having marginalised over the positions of the 
two clusters. 
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Figure 23. 1, 2 and 3<r \ 2 fit constraints on mi and mi for 2 clus- 
ters along line of sight, with noise appropriate for our space-based 
experiment, after including redshift measurements of cluster po- 
sitions. 

dius) at z = 0.6 behind another 8 x 10 13 M© cluster at 
z = 0.25, and applying a \ 2 fit for two masses and posi- 
tions using the method above. 

Figures 21 and 22 show the constraints we obtain on 
mass and position of the two clusters from the simulation, 
having marginalised over the 4-dimensional \ 2 distribution 
to find a subset of parameters. We see that, for this example, 
two configurations of clusters are possible: (a) two clusters 
very near each other in redshift; this is essentially the dis- 
covery of a 1-cluster solution, (b) One cluster at z = 0.25 
with another of similar mass at z = 0.6; this is the solution 
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Figure 24. Best fit model for <j> field along line of sight through 
centre of 2 clusters, m = 8 X 1O 14 M , at z = 0.25 and 0.5. The 
dotted line represents the <j> field arising from the z = 0.25 cluster 
alone. 



corresponding to our input scenario. In this case, we ob- 
tain greater accuracy in measuring the mass for the nearer 
cluster. This can be understood from Figure 24 (displaying 
masses 10 times larger for illustrative purposes): much of 
the background <j> amplitude is due to the first cluster, so 
the mass and position of this cluster is well-constrained; on 
the other hand, the second cluster's (j> contribution has less 
distance in which to rise, so this cluster's mass and position 
estimates are more affected by the noise amplitude. 

If this double solution were present when the fitting 
procedure was applied to real data, we could easily use our 
sample redshifts to confirm one scenario by looking for in- 
creased number-counts at the claimed cluster redshifts. By 
confirming the redshift of the cluster in this fashion, we can 
obtain better estimates of the mass, together with a more 
conclusive statement as to whether there are mass concen- 
trations in the background. This is demonstrated in Figure 
23. Here we have again simulated a cluster behind a cluster; 
we have now allowed the clusters to have known redshifts in 
our x 2 fit- In this case we find that a two-cluster fit is sub- 
stantially better than a one cluster fit (i.e. ni2 = 0) at the 2a 
level. We find best values for the masses of the clusters to be 
mi = (8.46lo 3g) x 1O 13 M and m 2 = (6.6311;^) x 10 13 M© 
(c.f. input mi, 2 = 8 x 10 13 Mq) within 2 arcmin radii. If we 
only constrain z\ , we find a best fit value for the background 
position of z-z = 0.6 ± 0.05 (c.f. input z = 0.6); thus if wc 
have found a low-redshift cluster, we can check if there are 
significant clusters behind it. 

We could increase the number of clusters in our x 2 fit, 
in order to seek for > 2 clusters along the line of sight. How- 
ever, in this case we will begin to fit the noise rather than real 
structures; this will be seen as too good a fit in our x 2 (i- e - 
A X 2 & y / 2(x/Ax)(y/Ay)(z/Az)). This restricts our ability 
to construct accurate 3-D maps using this method; however, 
we can at least map up to the first few significant mass con- 
centration in the z direction, and cosmological information 



could be gained by determining the statistical properties of 
the distance to the first or first few mass concentrations. 

One can envisage, therefore, a procedure consisting of 
(a) initial detection of mass concentrations using the 3-D 
distortion field alone, followed by (b) improved mass and 
background structure estimates by assigning accurate red- 
shifts from the visible matter associated with the detected 
foreground mass concentrations. 



8 CONCLUSIONS 

In this paper, we have developed and tested a practical 
method for 3-D reconstruction of the gravitational potential 
via weak lensing measurements, together with the more eas- 
ily obtained lensing potential. This methodology is based on 
the reconstruction equations of Kaiser & Squires (1993) and 
Taylor (2002), by which these local 3-D potentials can be 
calculated given knowledge of the lensing shear field in 3-D. 
This can be obtained by using shear estimators for galaxies 
with known redshifts. 

We have presented analytical forms for the shot-noise 
uncertainty in the convergence, lensing potential, gravita- 
tional potential and density field, noting that these fields 
become progressively more noisy for realistic survey sizes. 
We have also calculated the effects of only having finite sky- 
coverage for a survey and estimated the variance measured 
on such a survey, and the additional uncertainty in the re- 
construction due to structure beyond the survey boundary. 
In particular we have found that the contribution to the re- 
construction uncertainty in the differential lensing potential, 
A(j>, is dominated by large-scale structures. 

We have also shown that further sources of error, includ- 
ing photometric redshift errors, redshift-space distortions, 
and multiple scatterings of light rays, will not be dominant 
in our reconstruction process. 

In order to simulate the measurement of a 3-D gravita- 
tional field, we have calculated the expected lensing poten- 
tial due to a given mass field upon a 3-D grid. From this 
we have calculated the shear expected upon a galaxy im- 
age for a galaxy positioned anywhere in the 3-D grid. We 
have produced a catalogue of galaxies positioned in accor- 
dance with a realistic redshift distribution and given each 
an appropriate shear plus a random intrinsic ellipticity; this 
final galaxy shear catalogue was the information given to 
our reconstruction software. 

We have made reconstructions of the lensing and grav- 
itational fields, by smoothing the noisy shear data, calcu- 
lating the lensing potential according to Kaiser & Squires 
(1993) and using Taylor's equation (8) to find the gravita- 
tional potential. 

We have found that the method works well in recon- 
structing the full lensing and gravitational potentials in the 
absence of noise. However, the addition of Poisson sampling 
of the field at a finite set of galaxy positions, together with 
the noise due to the intrinsic ellipticity of objects, produces 
significant sources of error. We find that we obtain lensing 
potential maps with v ~ 6 in [3', 3', 0.05] bins in angle and 
redshift for a cluster of mass 8 x 1O 13 M . Unfortunately, cor- 
responding gravitational potential maps have only v ~ 0.5 in 
pixels of this size. However, applying Wiener filtering to this 
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gravitational potential, we can obtain u ~ 3 measurements 
of gravitational potential of clusters of mass 8 x 1O 13 M . 

This provides excellent prospects for obtaining cosmo- 
logically significant information directly from the measured 
gravitational potential field. For surveys with a redshift limit 
z — 1, mass concentrations J; 10 14 M Q can be directly 
mapped between 0.1 5^ z ^ 0.5, while statistical mapping 
can be used to examine the gravitational potential fluctua- 
tions on smaller mass scales. 

We have examined the effect of redshift uncertainties 
upon our simulations, and find that these contribute a much 
smaller error than the dominant intrinsic ellipticity noise 
term. 

Finally we have emphasised that, even on scales where 
a gravitational potential measurement is uncertain, the 3-D 
measurement of the <f> field is valuable, and can give us useful 
information about the statistics and topology of the mass 
field. In particular, we can obtain accurate measurements of 
the mass and position of clusters along the line of sight, and 
significantly detect the presence of clusters behind clusters. 

The methodology described here for engaging in 3-D 
gravitational mapping has numerous applications. We can 
obtain direct measurements of mass distributions in 3-D, 
which can act as an important cosmological probe via the 
mass function or cluster number counts. We can directly 
measure cross correlation functions between mass and light 
in 3 dimensions. Also, for possible 'dark' mass concentra- 
tions (e.g. Erben et al 2000), this 3-D mapping procedure 
allows us to measure both mass and radial position for such 
objects, which is quite impossible via conventional redshift 
methods. 



9 APPENDIX: CALCULATING 

UNCERTAINTIES ON THE 3D FIELDS 

Here we describe the details of our analysis for calculating 
uncertainties on the lensing and gravitational potentials. 



9.1 The lensing potential field 

In Section 3.1.2 we state the covariance of the 3-D lensing 
potential, 



(4>{r)4>{r'))sN — 49 2 d' 2 («(r) n(r' ')} s n ■ 



(66) 



We can calculate this covariance over the observed area, as- 
suming a flat sky, obtaining 

(<p(r)<f>(r'))sN = 

, 7r H , [ d 2 e"\n\e-e"\\n\e' -e"\6 D (r-r') (67) 
n 2 n{r)r 2 J A 

where the integral is taken over the survey area, A. The 
discrete case is an obvious change to a summation over bins. 
Transforming to the differential potential field, 



A0 : 



(68) 



where <f> is the mean field estimated from a finite survey of 
area A, 



>0(«), 



(69) 



we find the covariance (A0(r) Acj>(r')}sN is equivalent to 
equation (67) after transforming the kernel 

in \e - e'\ -> in |0 - 0'| - 1 / d 2 e in \e - e'\. (70) 

J A 

In the simple case of a circular survey with radius R we find 
i [ d 2 6'\n\0-0'\ = he/R) 2 -l+\nR 2 ]. (71) 

J A 

Note here we have assumed infinite resolution for the survey, 
or infinitely small pixels. 

The uncertainty in the 3-D lensing potential difference 
is then given by 



l/\A 2 {„\\ T™s 6 2 (0) 

(A</> (r))sN — 



■K 2 n(r) r 2 



8d(t — r'), 



where 

e 2 (0) 



If 



ln|0-0'l - 4- 



l2 6 ln|0-0'l 



(72) 



(73) 



which in general has to be evaluated numerically. In the 
special case of = and a circular aperture with angular 
radius 9, this can be evaluated analytically, giving 



(A0 2 (0 = 0))sn 



5 ~/L B 2 

24-7T n(r) r 2 



S D (r - r'), 



(74) 



where we have taken into account the conical geometry of 
the survey. This is the result discussed in Section 3.1.2. 



9.2 The Newtonian potential field 

In order to calculate the uncertainty on the 3-D Newto- 
nian potential, we must smooth the field; this is because 
we require a double differentiation of the lensing potential, 
but only sample this field at discrete points where there 
are galaxies. If we smooth in the radial direction with an 
arbitrary smoothing kernel, w(r), we find the resulting co- 
variance matrix of the Newtonian potential for a constant 
galaxy number density, n, in the distant observer approxi- 
mation, is 



7r 



e 2 (0,0') 



47r 2 n L 3 (r, r') 



where 
1 

L :i (r,r') ~ R 2 



r 2 r' 2 f 
R? J 



dyw"(r,y) w"(r',y) 



(75) 



(76) 



has units of inverse volume. Dashes on the window function 
denote derivatives with respect to distance. If the number 
density of sources is not a constant, these formulae must be 
altered by the substitution 



nL 3 (r, r') 



2/2 f 

r —^T j dyw"(r,y)w"(r', 



y)/n{y). 



(77) 



Again we can evaluate these expressions for the variance in 
a bin when r' = r and assuming a Gaussian weighting func- 
tion, w(r) = [\/27rY'||] _1 cxp(— r 2 /2r 2 ), with smoothing ra- 
dius r 1 1 . In this case to leading order, when r 3> r\ | , and tak- 
ing into account the conical geometry of the survey, the un- 
certainty on the Newtonian potential conveniently reduces 
to 
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2 5 7r ms / r \ (T\\ \ 2 2 Wittman D., Margoniner V. E., Tyson J. A., 2002, submitted to 



<$2(r))5iV " 6W^^P (^) (fl) ^ (?8) ApJL, astro-ph 0210120 

This is the result discussed in Section 3.1.3. 
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